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The Cherenkov Telescope Array (CTA) will be the next generation ground based observatory in 
very high energy gamma ray astronomy. The facility will achieve a wide energy coverage, start¬ 
ing from a threshold of a few tens of GeV up to hundreds of TeV by utilising several classes of 
telescopes, each optimised for different regions of the gamma-ray spectrum. The required energy 
resolution of better than 10-15% over most of the energy range and a goal of 5% systematic uncer¬ 
tainty on the measurement of the Cherenkov light intensity at the position of each telescope means 
that a very precise evaluation of the entire system will need to be made. The composite nature of 
the array means multiple camera technologies will be employed so multiple calibration systems 
and procedures will be necessary to meet the performance requirements. Additional constraints 
will come from the need to minimise observing time losses and that the observatory is foreseen 
to operate for tens of years, so both short and long term systematic changes in performance will 
need to be investigated and monitored. This contribution summarises the recommended camera 
calibration strategy of CTA based on the experience with current lACTs. 
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1. Introduction 

The camera of an imaging atmospheric Cherenkov telescope (lACT) consists of an array of 
pixels and their readout system. It begins with the photodetector that converts a Cherenkov photon 
to a photoelectron (p.e.) and amplifies this to a measurable current or voltage, which is in turn 
digitised and stored as digital counts (d.c.). The precision of the camera calibration will ultimately 
determine how accurately we can later reconvert the d.c. and use them to estimate the Cherenkov 
light intensity that arrived at the camera focal plane. Most of these calibration parameters will be 
characterised in the laboratory during the commissioning stage of the camera [1], but these initial 
values can vary due to ambient effects (e.g. temperature or night-sky background illumination 
level) or through ageing in the harsh operating environment. Here we refer to calibration as the 
on-site measurement and update of the lab. characterisation values, with measurements to be taken 
at intervals specific to the expected timescale of variation in the respective coefficient. 

The Cherenkov Telescope Array (CTA)^ [2] will be the next generation lACT facility, pro¬ 
viding the community with an open-access observatory for the observation of gamma rays with 
energies from a few tens of GeV to hundreds of TeV with unprecedented sensitivity and similar 
improvement in angular resolution and energy resolution. To achieve this large dynamic range in 
energy coverage the array will be comprised of multiple telescope sizes: a few 23 m diameter large 
size telescopes (LSTs) provide a low energy threshold; several 12 m diameter medium size tele¬ 
scopes (MSTs) provide an order of magnitude improvement in flux sensitivity in the 0.1 to 10 TeV 
region; and many 4 m diameter small sized telescopes (SSTs) will examine the lowest flux, highest 
energy events (E > 10 TeV). To meet the science case requirements the systematic uncertainty of 
the energy of a photon candidate (at energies above 50 GeV) must be < 15% and the energy reso¬ 
lution must be < 30% at 50 GeV and < 10% above 1 TeV; with the systematic uncertainty on the 
overall energy scale to not exceed 10%. The linear response of the CTA cameras must be recov¬ 
erable to < 5% up to 1000 p.e. per pixel. The MSTs and SSTs will also comprise of both single 
(-1M) and dual min'or (-2M) optics, leading to some very different camera configurations being 
used in the array. 

There are currently seven camera prototypes proposed for use in CTA, employing a wide 
variety of photodetector and data acquisition technologies: 

LST camera The LSTs will have AC coupled photomultipliers sampled by DRS4 Switched Ca¬ 
pacitor Arrays [3], this camera is optimised to lower the energy threshold by triggering a 
single telescope deep in to the night sky background (NSB) regime. A large dynamic range 
is achieved through the use of two gain channels. 

FlashCam A camera design for the MST [4] with DC coupled photomultipliers sampled by FADC 
and achieving a large dynamic range with non-linear gain pre-amplifiers of two distinct re¬ 
sponse regimes: a linear and a logarithmic region for small and large amplitudes respectively 
and a smooth transition between these regimes. 

NectarCam A camera design for the MST with AC coupled photomultipliers sampled by the 
NECTAr chip [5]. It also achieves its dynamic range thr'ough the use of two gain channels. 

* www.cta-observatory.org 
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SCT camera A camera for the dual-mirror Schwarzchild Couder Telescope (SCT) utilising Sili¬ 
con Photomultiplier (SiPM) readout by a TARGET ASIC [6] that samples the waveform [7]. 

SST-IM camera An adaptation of the FlashCam for the SST-IM, in this case DC coupled silicon 
photomultipliers (SiPM) are employed on the focal plane [8]. 

GCT camera A camera for a SST-2M. The photodetector is either Multi-Anode photomultiplier 
(MAPM) or SiPM that has the waveform sampled by a TARGET ASIC [9]. The camera 
saturates at about ~600p.e. but signal fitting (e.g. Time Over Threshold) can recover the 
signal to <5% bias up to ~1000p.e. 

ASTRI A camera for a SST-2M. AC coupled SiPM photosensors are read out by CITIROC which 
stores the peak of the electronically shaped signal in a pre-defined readout window width 
[10]. The dynamic range is achieved by having a low gain channel and a high gain channel. 

In light of this complexity it is necessary to have a clearly established set of guidelines for the 
calibration of the camera data to ensure that the CTA performance requirements are met, whilst si¬ 
multaneously minimising the complexity of the data analysis chain in order to simplify the software 
maintenance costs and make for easier upgrades. These guidelines will also assure the minimisa¬ 
tion in the complexity and fulfill the scientific requirements in the real time analysis [11] as well as 
offline in the data centre. 

2. Measurement of light intensity 

The response of a pixel (/) of solid angle D to a light of intensity I (photons per unit area and 
solid angle) is given by 

Vi{I) =kl + Pi = / • Aeff • Ruir • ^ • Tfj • rjQ i • rjcj ■ e ■ Gs,i • Geij + pt (2.1) 

where V (/) is the output in digital counts, possibly after some waveform processing and/or non¬ 
linearity corrections; Agff the effective mirror area after shadowing; /?Mir the focused mirror re¬ 
flectivity; Tf the transmission efficiency through the optical system, e.g. shadowing, lightcones, 
window, etc; rjQ the sensor quantum efficiency; rj^ the p.e. collection efficiency; e the electron 
charge; Gs the gain in the sensor; Gei the electronics conversion factor in terms of d.c. per p.e.; 
and p a baseline offset. For simplicity, we have omitted here the possible dependence of the factors 
on the angle of incidence, and assumed that quantities vary independently of the photon wave¬ 
length. In reality the dual mirror telescope designs will have photons striking the focal plane at 
much larger angles of incidence than the single mirror designs and also at different angles to any 
telescope mounted calibration light sources, which could introduce systematic differences of up to 
~ 30% if not properly accounted for. The degradation of the transmission and quantum efficiency 
coefficients is unlikely to be achromatic either and must again be carefully accounted for, espe¬ 
cially if there is any bias between the Cherenkov light spectrum from the air shower and that of the 
calibration light source (e.g. local muons). 

For CTA, the systematic error on the measurement of the absolute intensity (i.e. photons per 
square metre) of the Cherenkov light (post-calibration) at the position of each telescope must be 
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< 8%, with a goal of < 5%. The uncertainty in the estimate of the Cherenkov light intensity at 
each telescope is a combination of the uncertainties in the optical path efficiency, the Cherenkov 
photon detection efficiency af fhe focal plane and fhe sfafisfical and sysfemafic uncerfainfies in fhe 
phofon conversion, phofo-elecfron amplificalion and digitisation in fhe camera. Assuming a camera 
averaged response from eq. 2.1 and assuming fhere are no dominanf hofspofs fo bias fhe resulf we 
see fhaf fhe error {5S) in fhe Cherenkov image size {S, equivalenf fo fhe infensify) depends on 



where 8V is fhe uncerfainfy in fhe electronics signal, which contains the photon conversion ef¬ 
ficiency, and uncerfainfies in knowledge of fhe precise amplification and digifizafion; 8p is fhe 
uncerfainfy from fhe pixel baseline; 8FF is fhe uncerfainfy in fhe camera flaf-fielding; 8T is fhe 
uncerfainfy in fhe lighf fransmission; and 8MC is fhe uncerfainfy infroduced by approximafions 
used in fhe Monfe Carlo simulafion. If we divide fhe error budgef equally befween fhe phofon 
fransmission before reaching fhe camera (i.e. mirrors, lighfguides, window and shadowing) and 
from fhe conversion fo p.e. onwards (i.e. phofodefecfion, amplificafion and digifisafion) we can 
reach fhe goal requiremenf allowing 1% uncertainly in fhe digifisafion, 0.2 p.e. uncerlainly in fhe 
pixel baseline measuremenf and 2% uncerlainly in fhe flal-fielding of fhe camera response. Some, 
like fhe uncerlainly in fhe pedeslal baseline, have already been shown fo be achievable on cur- 
renl generalion inslrumenls; olhers, like measuring fhe uncertainly on fhe fransmission, are a slep 
furlher in fhe calibralion demands compared fo currenf inslrumenls (due fo fhe addifional need fo 
conlrol fhe wavelenglh and angular dependence aspecls); olhers slill, like flal-fielding fhe curved 
focal plane of fhe dual mirror lelescope syslems from shorl illuminalion disfances have presenfed 
challenges never before faced in Ihis field and so novel techniques musl be Ihoughf of (e.g. [12,13]). 
If fhe chromalic biases befween Cherenkov lighf from muons and from air showers can be kepi fo 
a minimum fhen an end-fo-end calibralion using muons rings wilh a precision of < 5% can be mef, 
buf resolving fhe individual confribulions can still be challenging. 

Tracking where issues in fhe sysfemafic uncertainties may arise is best achieved when multi¬ 
ple, or redundant, methods are at ones disposal. Figure 1 shows a mock data set for the relative 
efficiency of a generic camera over a 5 year period as measured by the multiple methods of calcu¬ 
lating the p.e./d.c. value from single photoelectron runs and from using muon rings, and the relative 
measurements of cosmic-ray throughput [14] based on image size and from triggering pixels only. 
Each method is sensitive to different levels of systematic effects that change the expected number 
of digital counts received from the camera and in turn to how the system is simulated to produce 
lookup tables for data analysis. To exaggerate the effects for display purposes the camera is sim¬ 
ulated with a 20% drop in pixel efficiency per year through ageing (compensated for by raising 
the voltages once per year); a 25% change in the Cherenkov light yield due to seasonal variations 
in atmospheric density and a 10% drop in transmission due to an enhanced high altitude aerosol 
layer (such as from calima or biomass burning) in month 6 and 7 of a year; and a 15% change 
in the light focused onto a pixel occuring 20 months after the start of observatory operations, e.g. 
as might be expected if the telescope were refocused between infinity and the altitude of shower 
maximum. Only at very specific times do all the methods give consistent results. The single pho¬ 
toelectron measurement will only be sensitive to changes in the camera specific hardware, whereas 
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the other end-to-end measurements are also susceptible to the systematics of the atmosphere and 
the telescope structure. The seasonal changes in the atmosphere at certain times can counteract the 
camera hardware ageing, at others enhance them. The change in focus affects measurements that 
use the per pixel amount of light, such as muon rings and triggering pixels, but not ones that use all 
of the image light, hence the splitting between the different throughput measurements at month 20. 
Where achievable we recommend a primary method for the calculation of calibration parameters 
and secondary (or tertiary) methods that can be used as cross-checks on a less frequent basis, or in 
the event of failure of hardware for the primary method. An additional consideration to take into 
account is that a minimal amount of time must be taken from the observing schedule to deal with 
calibration issues in dark time. As such, any calibration that can be performed simultaneously with 
the acquisition of observational data must do so, this will also have the added benefit of giving 
calibration information that is relevant for the current ambient conditions at the time of data taking. 



Figure 1: Mock up of a telescope’s relative performance over a 5 year period as measured by four different 
methods: single p.e. (black cross); cosmic-ray throughput measured by image size (red star); cosmic- 
ray throughput measured by 2nd brightest pixel (green circle); muon rings (blue cross). Degradation is 
exaggerated to clearly see the differences in the figure. The error bars are for a 2.5% random statistical error 
that has been included for the data points. Systematic effects including PMT ageing (red line), seasonal 
variation in Cherenkov light yield (green line) with an additional aerosol layer in months 6 & 7 of a year and 
periodic changes including a yearly increase in the gain and a refocusing of the camera at month 20 are also 
plotted. 


3. Calibration hardware and procedures 

Hardware, procedures, methods and interfaces for the following calibration must be in place, 
though not every one will be relevant to all camera designs. They fall in to the categories of regular 
monitoring (> once per night, e.g. [15]) and long term monitoring (at least once per year, that 
check the coefficients measured from the camera characterisation stage [1]). 

3.1 Regular Monitoring: on telescope systems 

Pixel baseline and its fluctuations Should be determined under actual observing conditions in or¬ 
der to account for potential baseline shifts due to background illumination (depending e.g. 
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on whether electronics is AC or DC coupled) and ambient environmental effects (e.g. tem¬ 
perature). Various well-proven techniques exist to determine pixel baseline, e.g.: random- 
triggered events can be injected into the data stream (this is the preferred scheme for CTA); 
alternatively, baselines can be determined from blank regions of shower images; or in the 
case of waveform readout, baselines can be determined from samples preceding the trig¬ 
gered signal, provided a sufficiently wide readout window is used. The electronics baseline 
(i.e. independent of NSB contribution) can be determined by force triggering the camera 
with the camera lid closed. 

Relative pixel response to incident light (intensity flat fielding) Flat-fielding of fhe camera can 
be achieved by fwo mefhods: illumination of fhe camera wifh a lighf pulser wifh a lighf fronl 
of known intensify variation (preferably uniform); or fhe pixel response for muon rings. The 
lighf pulser will usually be fhe simples! and mosf effecfive solution. 

Relative pixel response to incident light (timing flat fielding) The timing offsets between pixels 
should be identified to sufficient precision using the intensity flat-fielding events. 

Effective number of photoelectrons / excess noise Several techniques can be employed to de¬ 
termine the gain of a photosensor: the fluctuation in the response to large calibration light 
pulses (see e.g. [16, 17]), which is the easiest scheme to implement on a nightly basis whilst 
minimising the impact to the observation schedule by using the flat-fielding events; exposure 
of the camera pulsed light at low level, under suitably dark conditions (i.e. using a shel¬ 
ter or camera mask), such that the single photoelectron signal can be directly determined, 
which can be used as an accurate cross-check on a less frequent basis, or comes for free 
with SiPM based cameras; exposure with night-sky background, using timing information to 
select pixels which have a signal likely related to one or few night-sky photoelectrons. 

On telescope light flasher systems can be simple, and serve primarily for flat-fielding. LED 
systems with pulse width and spectrum roughly matched to the Cherenkov spectrum seem most 
promising, following e.g. the concept of [12], though the exact lACT requirements may still require 
the use of a laser based system, e.g. [18, 4]. They should provide a known illumination light front 
shape of the camera with a fixed pulse height of > 30p.e. and preferably 100 p.e.. Strict long-term 
stability of the flash amplitude is not required, short-term stability (minute-scale) must be good 
compared to the statistical variation of the pixel signals. It should be possible to trigger the flasher 
during a run to inject calibration events without significantly affecting observations. 

3.2 Long Term Monitoring: high-performance light flasher systems 

A high-performance light flasher system should be available to characterize telescopes at 
longer intervals (e.g. once per year). The flasher system can either be placed on the lACT, or 
- possibly simpler but requiring larger light intensity - a movable system that is placed at a set 
distance from the telescope (e.g. [13]). In the latter case, it should reasonably uniformly illuminate 
the whole mirror surface. The system should ideally provide the following features: 

• Combination of pulsed and DC light 
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• DC light with selectable wavelength over the range of sensor sensitivity, and adjustable in¬ 
tensity up to roughly full-moon equivalent. Absolute calibration of the DC light intensity is 
not required. 

• Pulsed light with selectable wavelength over the range of sensor sensitivity, and adjustable 
intensity up to saturation of the sensors or electronics chain. Pulse length should match 
Cherenkov pulse length. Relative light level should be calibrated with a precision of better 
than 3%, e.g. using a precision optical attenuator. Absolute light level should be calibrated 
with a precision of 5%, e.g. using calibrated photodiodes at some distance. 

• Built-in battery power and wireless control would increase flexibility. 

The system will then be used to measure 

Pixel linearity To determine pixel response at several intensity levels up to and beyond 1000 p.e., 
particularly for dual gain channel or non-linear amplification systems. 

Pixel spectral response Spectral response needs to be measured in the lab and occasionally be 
verified in the field. Spectral response must be known sufficiently well that the calibration 
reference light source yield relative to the air shower light yield and night-sky background 
response is known with a precision at a level of ~ 1 — 2% - this is especially true for calibra¬ 
tion methods using the Cherenkov light emitted by muons. A calibration system providing 
four colours ranging from UV to red should be at least available once per year. 

Pixel response for different light incidence angles An appropriate solution would be to have a 
high-performance angular response calibration system, consisting of a single pe illumination 
source that can illuminate at a range of angles from face-on to the maximum expected ray 
angle (up to 70° for the dual mirror telescope systems), available in the on-site lab to calibrate 
a representative number of pixels off-telescope, and a system placing a remote light source 
at about twice the focal distance, with the telescope parked and the reference light front thus 
following a similar optical path via the mirrors to the focal plane; however, shadowing issues 
and reproducibility of the method must be demonstrated beforehand. 

Integration Window Special calibration runs with the extraction window at maximum can be 
compared to runs with the standard shorter operating window. Flat-fielding events with a 
similar FWHM to Cherenkov events can be used to determine an appropriate window size 
relatively quickly, but Cherenkov images will have a larger time structure across the image 
which can result in signal loss. It is recommended to determine these losses through specific 
runs with the telescopes pointing at zenith and comparing the image size distributions as a 
function of the reconstructed impact distance. 

4. Summary 

The ambitious science requirements and long operational life span of CTA, combined with the 

increased complexity of operating multiple telescope types, means that clearly defined calibration 

procedures in the lACT camera calibration need to be in place if systematic uncertainties are to be 
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kept to a minimum. The camera calibration can be split according to whether a coefficient needs 
to be regularly or infrequently measured: the former preferably done with events that are taken in 
conjunction with science observations; the latter that can be done with more specialised hardware 
in dedicated campaigns. 

We gratefully acknowledge support from the agencies and organizations listed under Funding 
Agencies at this website: http://www.cta-observatory.org/. 
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